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The  report  describes  three  position  estimation  procedures.  The 
first  is  for  estimates  based  on  bearings  taken  on  or  from 
stations  with  uncertain  locations.  The  second  is  for  use  with 
two  or  more  lines  of  position.  The  third  is  for  combining 
estimates  from  different  sources.  The  initial  version  of  the 
report  contained  an  error  in  a  relation  relating  to  the  third 
I  procedure.  Appendicies  to  the  report  contain  a  description 
and  listing  for  a  program  that  can  be  used  to  implement  the  first 
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This  revised  report  contains  a  correction  to  the  expression  for 
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at  the  bottom  of  page  14.  In  addition,  it  contains  a  development 
in  an  added  Appendix  1  of  a  procedure  associated  with  the  position 
estimation  procedures  discussed  in  Sections  II  and  III  of  the 
report.  It  also  contains  a  description  and  a  program  listing  in 
an  added  Appendix  2  for  a  program  called  PEST  that  can  be  used  to 
implement  the  first  two  position  estimation  procedures  and  for  a 
program  called  COMP  that  can  be  used  to  implement  the  position 
estimation  procedure  discussed  in  Section  IV  of  the  report.  Both 
programs  are  written  in  BASIC  for  a  Radio  Shack  TRS-80  Model 
100  Portable  Computer. 


I.  Introduction 


This  report  describes  position  estimation  procedures  that 
are  based  on  models  which  relate  positional  uncertainty  to 
various  measurement  and  estimation  errors.  The  procedures  were 
developed  to  be  used  in  analyses  that  relate  the  effect  of  posi¬ 
tional  uncertainty  on  tactical  performance  through  such  factors 
as  weapon  accuracy  as  well  as  to  be  used  operationally. 

In  the  models,  positions  are  on  a  plane  surface  (flat  earth 
model) .  Because  of  this  condition,  the  models  are  not  intended 
for  use  in  situations  where  the  earth's  figure  is  significant. 

In  addition,  positional  errors  are  determined  by  independent 
normally  distributed  random  variables  with  known  means,  variances 
and  covariances.  The  support  for  this  condition,  other  than  its 
mathematical  convenience,  is  that  it  has  been  used  by  others, 
for  example,  see  Reference  1.  To  use  the  models,  one  is  required 
to  specify  the  means,  variances  and  covariances  of  the  error 
random  variables.  For  these  models,  this  can  be  done  by  speci¬ 
fying  the  systematic  errors  (biases)  and  the  error  magnitudes 
(standard  deviations) . 

The  procedure  that  is  described  in  Section  II  relates  position 
estimates  that  are  based  on  bearings  on  or  from  stations  to 
station  position  uncertainty.  One  application  for  the  model  is 
the  analysis  of  the  effect  of  sonobuoy  position  uncertainty  on 
position  estimates  determined  with  passive  directional  sonobuoys. 

The  procedure  that  is  described  in  Section  III  relates  posi¬ 
tion  estimates  that  are  based  on  lines  of  position  to  line  of 


position  uncertainty.  In  addition  to  estimating  the  uncertainty 
in  position  estimates  based  on  celestial  observations,  the  pro¬ 
cedure  could  be  used  to  determine  error  ellipses  for  LORAN  fixes 
if  the  standard  deviation  values  required  by  the  procedure  could 
be  obtained. 

The  procedure  that  is  described  in  Section  IV  can  be  used 
to  combine  position  estimates  from  various  sources.  The  proce¬ 
dure  v’hich  is  based  on  conditions  that  should  not  be  too  restric 
tive  in  most  cases  provides  both  a  composite  position  estimate 
and  error  ellipse. 


II.  station  Position  Uncertainty  and  Position  Estimates 

The  procedure  that  is  described  in  this  section  relates 
position  estimates  that  are  based  on  bearings  on  or  from  stations 
to  station  position  uncertainty.  The  procedure  is  based  on  a 
model  that  is  an  extension  of  one  that  is  described  in  Appendix 
1.  The  model  is  defined  as  follows:  Each  station  position  error 
is  determined  by  an  independent  bivariate  normal  distribution 
with  a  zero  mean  vector  and  a  known  covariance  matrix.  Observed 
bearing  lines  on  or  from  a  station  are  parallel  to  true  bearing 
lines.  The  distance  of  each  observed  bearing  line  from  its  true 
bearing  line  is  determined  by  an  independent  normal  distribution 
witu  a  zero  mean  and  a  standard  deviation  o. 

Becuuse  a  station's  position  error  is  determined  by  a  bivar¬ 
iate  normal  distribution,  the  perpendicular  distance  s  between  a 
line  at  the  assumed  location  of  an  observed  bearing  line  and  the 
observed  bearing  line  is  determined  by  a  normal  distribution. 

The  relation  between  the  bivariate  normal  distribution  that 
describes  the  station's  position  uncertainty  is  indicated  in 
Figure  1.  In  the  figure,  the  positive  y-axis  direction  is  north, 
the  positive  x-axis  direction  is  east,  and  the  origin  of  the 
coordinate  system  is  at  the  assumed  station  position.  The  x'y'- 
coordinate  system  is  oriented  so  that  the  positive  x'axis  is 
coincident  with  the  major  axes  of  the  elliptical  contours  of 
the  bivariate  normal  distribution  that  determines  the  station 
position  error  and  so  that  the  bearing  6  of  the  positive  x'-axis 


satisfies  the  condition:  0°  <  6  <  180°.  The  angle  6  is  an 
observed  bearing.  The  x"y "-coordinate  system  is  oriented  so 
that  the  positive  y"-axis  is  in  the  direction  of  the  observed 
bearing  line.  As  a  consequence  of  these  relationships, 

2  2  2  2  2 
0^  =0^,  sin  (6-0)  +  o„,  cos  (6-0) 

s  X  y 

2  2 

where  and  are  the  elements  of  the  station  position 

error  covariance  matrix  relative  to  the  x ' y ’ -coordinate  system. 

A  procedure  for  determining  position  estimates  that  are 
based  on  bearings  on  or  from  stations  at  known  positions  is 
described  in  Appendix  1  of  this  report.  The  procedure  is  based 
on  a  model  that  is  equivalent  to  one  that  is  described  in  Ref¬ 
erence  2.  The  procedure  that  is  described  in  this  section  is 
based  on  a  model  that  is  an  extension  of  it. 

The  procedure  in  Appendix  1  of  this  report  is  based  on  a 
model  in  which  a  station's  bearing  error  is  determined  by  a 
normal  distribution  with  zero  mean  (bias)  and  standard  devia¬ 
tion  e.  The  bearing  error  is  related  to  the  distance  on  a  cir¬ 
cular  arc  between  a  station's  true  bearing  line  and  observed 
bearing  line.  The  arc  is  on  the  circle  with  its  center  at  the 
station  that  passes  through  an  initial  estimate  of  an  object's 
position.  This  distance  is  determined  by  a  normal  distribution 
with  mean  zero  and  standard  deviation  o  =  re  where  r  is  the 
range  of  the  initial  estimate  from  the  station  and  the  standard 
deviation  e  is  measured  in  radians.  In  the  model,  arc  distance 
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is  approximated  using  a  first  order  approximation  which  in 
effect  replaces  the  circle  with  its  tangent  line  at  the  initial 
estimate.  As  a  consequence,  the  distance  u  on  the  tangent  line 
between  the  observed  bearing  line  and  the  true  bearing  line  is 
determined  by  a  normal  random  variable  with  mean  zero  and  stan¬ 
dard  deviation  o.  The  distance  u  can  be  expressed  in  terms  of 
w,  the  distance  on  the  tangent  line  between  the  observed  bearing 
line  and  the  initial  estimate  that  is  also  determined  by  a  normal 
random  variable  with  standard  deviation  a,  and  v,  the  distance 
on  the  tangent  line  between  the  true  bearing  line  and  the  initial 
estimate.  And,  as  shown  in  Appendix  1  of  this  report,  this 
distance  can  be  expressed  in  terms  of  the  unknown  coordinates 
of  the  object's  position. 

The  effect  of  station  position  uncertainty  is  accounted  for 
by  the  distance  s  between  the  observed  bearing  line  and  the 
assumed  observed  bearing  line.  In  the  model  s  is  determined  by 
a  normal  distribution  with  mean  zero  and  standard  deviation  a 

s 

as  given  above.  This  approximation  is  consistent  with  the  first 
order  approximation  of  arc  distance.  As  a  consequence  of  these 
two  approximations,  all  of  the  bearing  lines  are  replaced  by 
lines  parallel  to  the  line  joining  the  initial  estimate's 
position  and  the  station's  assumed  position,  both  of  which  are 

known  positions.  The  geometry  involved  in  shown  in  Figure  2. 

I 

The  modified  relationships  resulting  from  the  introduction  of 
station  position  uncertainty  are  shown  in  Figure  3.  The 

I 
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assumed  location 
of  the  observed  bearing  line 


The  quantities  u,  v  and  w  correspond  to  the  quantities 
u,  V  and  w  in  Appendix  1  of  Reference  3.  The  auxiliary 
quantities  u'  and  v?'  are  defined  in  the  figure.  The 
replacement  of  by  +  Og^  in  the  procedure  in 
Appendix  1  of  this  report  is  justified  by  noting  that 
with  this  substitution  u'  and  w*  are  equivalent  to  u 
and  w. 


important  thing  to  note  is  that  u*  =  u  -  s  is  determined  by  a 

normal  distribution  with  mean  zero  and  standard  deviation 
2  2  h 

[c  +  Og  ]  and  that  otherwise  U'  is  equivalent  to  U  with 

respect  to  the  procedure  in  Appendix  1  of  this  report.  As  a 

consequence  of  this,  the  procedure  can  be  extended  to  include 

2  2-5 

station  uncertainty  by  replacing  a  by  [a  +  Cg  ] ^  where  ever 
it  is  used.  In  this  case,  o  =  re  where  r  is  the  range  of  an 
initial  estimate  of  an  object's  position  from  the  assumed  station 
position  and  e  is  the  bearing  error  (standard  deviation)  in 
radians  of  the  bearings  associated  with  the  station. 


III.  Line  of  Position  Uncertainty  and  Position  Estimates 

The  procedure  that  is  described  in  this  section  relates 
position  estimates  that  are  based  on  lines  of  position  to  line 
of  position  uncertainty.  The  procedure  is  based  on  a  model  that 
is  defined  as  follows:  Lines  of  position  are  straight  lines. 
Observed  lines  of  position  are  parallel  to  true  lines  of  posi¬ 
tion.  The  distance  of  an  observed  line  of  position  from  a  true 
line  of  position  is  determined  by  an  independent  normally  dis¬ 
tributed  random  variable  with  known  mean  and  standard  deviation 
Lines  of  position  are  specified  in  terms  of  a  rectangular  coor¬ 
dinate  system  with  the  origin  at  a  reference  point  as  shown  in 
Figure  4.  For  celestial  navigation,  an  appropriate  choice  for 
the  reference  point  would  be  the  assumed  position.  Since  bear¬ 
ing  lines  are  lines  of  position,  this  model  differs  from  the 
model  that  is  described  in  Section  II  only  in  terminology. 
However,  operationally  the  use  of  the  model  that  is  described 
in  this  section  differs  in  the  way  the  standard  deviation  of  the 
distance  of  the  line  from  a  true  line  is  determined.  The  stan¬ 
dard  deviation  associated  with  each  line  of  position  must  be 
specified.  If  this  can  be  done,  the  procedure  can  be  used. 

As  an  example,  suppose  that  the  values  are  for  the  first 

line  of  position  and  02  for  the  second  where  >  O2  and 
that  the  lines  of  position  intersect  at  a  90°  angle.  In  this 
example,  the  minimum  area  confidence  (probability)  region  is 
an  ellipse  that  is  centered  on  the  estimated  position.  And  the 
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Figure  4.  The  geometry  when  only  two  lines  of  position  are  used 
with  the  procedure  in  Appendix  1  of  this  report.  The  lines  of 
position  correspond  to  bearing  lines  for  observed  bearings  and 
©2  with  respect  to  the  procedure.  The  procedure  determines  the 
size  and  the  orientation  of  elliptical  confidence  regions  of  con¬ 
fidence  p,  the  location  of  the  initial  estimate  and  the  location 
of  the  final  estimate.  The  locations  of  the  estimates  correspond 
when  only  two  lines  are  used.  If  the  lines  were  obtained  from 
sextant  observations,  the  assumed  position  should  be  the  reference 
position  in  which  case  the  lines  of  position  would  be  determined 
by  the  azimuth  anges  a.  and  and  the  distances  p,  and  p_. 


estimated  position  is  at  the  intersection  of  the  bearing  lines. 

2 

For  a  confidence  (probability)  of  containment  of  1  -  exp(-k  /2) , 
the  major  axis  of  the  ellipse  is  coincident  with  the  first  bear¬ 
ing  line  and  it  is  of  length  2kOj^,  the  minor  axis  of  the  ellipse 

is  coincident  with  the  second  bearing  line  and  it  is  of  length 

2 

2ko-  and  the  area  of  the  ellipse  is  Trk  . 


IV.  A  Composite  Position  Estimate 

The  procedure  that  is  described  in  this  section  is  for 
combining  position  estimates  for  an  object  from  independent 
sources.  It  is  based  on  the  following  model:  The  rectangular 
coordinates  of  each  position  estimate  are  determined  by  an  inde¬ 
pendent  bivariate  normal  distribution  whose  covariance  matrix 
is  known  but  whose  mean  vector  is  not  known.  The  components  of 
the  mean  vector  for  each  of  the  distributions  are  the  unknown 
coordinates  x  and  y  of  the  object.  This  model  implies  that  the 
natural  logarithm  of  the  likelihood  function  for  a  set  of  n 
estimates  can  be  expressed  as  follows: 

n  ^  ,  -1^ 

log  L  =  K  -  1/2  Z  (x^  -  x)  (Xj^  -  x) 
i=l 

where  K  is  a  constant,  x^  is  an  estimate  vector  with  components 
and  x  is  the  common  mean  vector  with  components  x  and  y, 

the  unknown  coordinates  of  the  object,  and  E^  is  the  covariance 

2  2  ... 
matrix  with  elements  a.^  ,  a . and  .  The  maximum  likeli- 

IX  ixy 

hood  estimates  x  and  y  of  the  unknown  coordinates  x  and  y  are 
the  solutions  of  the  two  simultaneous  linear  equations  deter¬ 
mined  by 


9 ( log  L) 


9 ( log  L) 
9y 


The  equations  can  be  written  as 


A  X  +  B  y  =  D 
B  X  +  C  y  =  E 


and  their  solution  as 


X  = 


CD  -  BE 


-  AE  -  BD 

y  =  - 71 


AC  -  B  AC  -  B 

where  A  =  Za^^,  B  =  Ib^,  C  =  Zc^,  D  =  Z{aj^x^  +  b^y^)  , 


E  =  +  c.y.),  a.  .  o^-^/d.,  =  -  OiWd^,  c.  = 


2  2  2 

d-  =  0.^  o.-  -  0.;'^  and  all  of  the  sums  are  for  i  from  1  to  n 
1  IX  ixy 


since  they  are  linear  combinations  of  the  estimates  x^  and 

y^,  the  estimates  x  and  y  are  determined  by  a  bivariate  normal 

distribution.  Consequently,  all  that  is  required  to  determine 

this  distribution  is  its  mean  vector  with  components  and  y 

2  2 

and  its  covariance  matrix  with  elements  and  The 

X  y  xy 

mean  vector  is  determined  by 


y-  =  f{(CD  -  BE) /(AC  -  B^)}  =  X 

and 

y-  =  E{(AE  -  BD)/(AC  -  B^)}  =  y 


And,  the  covariance  matrix  is  determined  by 

=  E  { (CD  -  BE)  -  E  (CD  -  BE) }^/(AC  -  B^)^ 

=  {C^ (F+I  +  2L)  -  2CB(H+K+M)  +  B^ (G+J  +  2N)  }/ (AC  -  B^)^ 


<>i 


a-  =  t{(AE  -  BD)  -  fc(AE  -  BD)}  /(AC  -  B  ) 

=  {A^(G+J+2N)  -  2AB(H+K+M)  +  B^ (F+I+2L) }/ (AC  -  B^)^ 


^xy  ^  E{[ (CD-BE) -E (CD-BE) ]• [ (AE-BD) -E (AE-BD) ] }/(AC  -  B  ) 

=  { (AC+B^) (H+K+M) -CB(F+I+2L)-BA(G+J+2N) }/(AC  -  B^)^ 

where  F  =  ,  G  =  H  =  Ea.b.o.^^  ,  I  =  Eb.^o.^^ 

1  IX  '  1  IX  1  1  IX  '  1  ly 

2  2  2  2 
J  =  Ec.  a.-'  ,  K  =  Eb.c.a.-'  ,  L  =  Ea.b.o.-^,  M  =  E(a.c.+b. 

1  ly  '  1  1  ly  '  11  ixy'  i  i  i  ixy 

and  N  =  Eb.c.a.^''  where  all  of  the  sums  are  for  i  =  1  to  n. 

1  1  ixy 

By  using  arguments  given  in  Appendix  1  of  this  report,  one 

can  show  that  the  axes  of  the  elliptical  confidence  regions 

associated  with  x  and  y  are  coincident  with  an  x'y '-coordinate 

system  where  the  transformation  from  the  xy-coordinate  system 

to  this  system  is  the  coordinate  axes  rotation  through  the 

2  2 

angle  >  defined  by  tan  2y  =  2o''''/(o^  -  o-'  )  .  For  a  confidence 

xy  y  X 

p,  from  Reference  3,  the  minimum  area  confidence  region  is  an 

2 

ellipse  with  semi-axes  k  and  k  o(',.  ,  and  area  irk  ooiCJ('i 

X  y  /  ^  y 

where  k  =  [-2  In(l-p)]*^  ,  > 

2  2  2  2  2 

2  2  2  2  2 

o^,  =  0''  sin  Y  +  20''-'  cos  y  sin  y  +  o;;  cos  y 
y '  X  xy  y 

and  the  center  of  the  ellipse  is  at  the  point  (x,y) .  In  this 
coordinate  system,  ocir. i  =  0  • 


The  above  equations  can  be  used  to  specify  a  composite  posi 

tion  estimate  in  terms  of  the  location,  orientation  and  size  of 

an  elliptical  confidence  region  which  is  generally  the  form  in 

which  position  estimates  of  the  kind  that  are  being  considered 

2  2 

here  are  specified.  But,  since  values  of  and 

X  y  xy 

are  required  for  each  of  the  n  estimates  that  is  being  combined 
a  way  is  needed  for  determining  these  values  given  the  orienta¬ 
tion  and  size  of  an  elliptical  confidence  region.  A  procedure 
to  do  this  when  the  orientation  is  given  in  terms  of  the  direc¬ 
tion  6  of  the  major  axis  and  the  size  is  given  in  terms  of  the 
lengths  SMJ  and  SMI  of  the  semi-major  and  semi-minor  axes  and 
the  confidence  p  is  described  next. 

By  using  an  xy-coordinate  system  in  which  the  positive 
y-axis  direction  is  north  and  the  positive  x-axis  direction  is 
east  and  with  the  convention  0°  <  6  <  180®,  the  dependence  of 
the  value  of  the  rotation  angle  Y  and  of  the  order  relation 

between  o^,  and  ,  on  the  value  of  the  major  axis  direction  6 
X  y 

is  indicated  by  the  following  table: 

0°<6<  45°:  Y  =  <5  and 

—  y  X 

45°  5  6  <  135°:  Y  =  6  -  90°  and  o^,  >  o^, 

135°  <  6  <  180°:  y  =  <5  -  180°  and  o-,  >  o-, 

With  an  order  relation  and  a  value  for  p,  values  for  ,  and 

Oy,  can  be  determined  with  values  for  SMJ  and  SMI.  With  values 

2  2 

for  »  cj-',  and  y»  values  for  o-  ,  and 

X  y  X  y  xy 

from  the  following  equations: 


can  be  found 


X  ~  x'  y'  ' 

2  2  2  2  2 
a''  =  O'-,  sin  y  +  o^i  cos  y 

y  x'  y 


2  2 

=  (O'',  -  O'-  )  sin  Y  cos  Y 

xy  y '  X ' 

which  can  be  obtained  by  inverting  the  equations  above  for  a^, 
2 

and  o^, 

As  an  example,  suppose  the  data  in  the  following  table 
represent  three  independent  position  estimates; 


X 

y 

6 

SMJ 

SMI 

-3.7 

18.1 

59“ 

36 

20 

11.8 

00 

• 

105“ 

37 

11 

0 

0 

146“ 

45 

23 

Here,  distances  are  in  nautical  miles  and  p  =  .86  in  each  case 
For  this  example  with  values  in  square  nautical  miles: 

®lx^  =  264.58,  =  159.42  and  =  98.89 

=  321.35,  =  51.15  and  =  -78 

2y  2y  2xy 

^3x^  =  249.20,  =  389.30  and  =  -173.4 


These  values  give  the  following  composite  estimate: 

X  =  -2.46  nautical  miles  and  y  =  12.26  nautical  miles. 
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6.21  nautical 


For  this  case,  o-  =  10.28  nautical  miles,  = 

miles  and  =  -55.1  square  nautical  miles.  And,  for  k  =  2; 
xy 

SMJ  =  20.56  nautical  miles,  SMI  =  12.42  nautical  miles  and 
6  =  119°.  The  composite  confidence  region  and  its  three  com¬ 
ponent  confidence  regions  are  shown  in  Figure  5. 

As  a  second  example,  suppose  each  position  estimate  is  deter 

mined  by  a  circular  normal  distribution.  Then  o.^  =  o.^  =  a. 

XX  X  y  X 

and  =  0  for  i  =  1  to  n.  In  this  case,  the  composite 

estimate  is: 

X  =  / (Ll/a^)  ,  y  =  (Zy^/o^) / (El/o^)  , 

=  n/(Zl/a^)^  ,  =  n/(Zl/a^)^  and  =  0. 

In  this  example,  since  x  and  y  are  determined  by  a  circular 
normal  distribution,  the  minimum  area  confidence  regions  are 
circles  and  orientation  is  not  an  issue. 


Figure  5.  The  ellipses  define  the  confidence  regions  of 
the  first  example.  The  composite  confidence  region  is  in 
black.  The  position  estimates  are  at  the  center  of  the 
ellipses.  The  numbers  indicate  the  order  of  the  estimate 
in  the  table  on  Page  17. 


APPENDIX  1: 


A  Position  Estimation  Procedure 


Consider  a  target  whose  position  is  unknown  and  a  set  of 
reference  stations  whose  positions  are  known.  Assume  conditions 
are  such  that  observed  target  bearings  from  the  stations  can  be 
considered  to  be  values  of  independent  normal  random  variables 
whose  means  are  equal  to  the  true  target  bearings  and  unknown  but 
whose  standard  deviations  are  known.  The  position  estimation 
procedure  that  is  described  here  is  a  maximum  likelihood  estimation 
procedure  that  is  based  on  this  assumption  and  the  assumption  that 
the  conditions  are  such  that  the  stations  and  the  target  can  be 
considered  to  be  located  on  a  plane  tangent  to  the  earth's 
surface  at  a  point  in  their  neighborhood. 

With  the  above  assumptions,  the  likelihood  of  observed 
bearings  0^,  9^/  •••,  9j^  from  stations  labeled  1,  2,  n  is: 


0  )  = 
n 


exp  - 


>2.2 
> .  /e . 

1  1 


where  •••>  'I’pj  tdie  unknown  station  true  bearings  and 

e,,  e-,,  ...,  e  are  the  known  stations  standard  deviations. 

1  2  n 

To  a  first  order  approximation,  the  set  of  bearing  estimates 
;  that  are  determined  by  the  procedure  make 

L(i3^,  ...,  0^)  a  maximum  subject  to  the  constraint  that  the 

bearing  lines  determined  by  a  set  of  bearing  estimates  must  all 
pass  through  a  common  point.  The  common  point  is  the  estimate 
of  the  target's  position. 


2) 


In  order  to  impose  the  constraint  on  the  estimate  through  4'^/ 

consider  the  quantities  u.  =  r.(0.  -  4>  ■  )  ,  v.  =  r.(4>-  “  3.)  and 

-  B^)  where  i  represents  any  station  number  from  1  to  n- 

Here,  is  the  target's  true  bearing,  is  the  bearing  of  an 

initial  estimate  of  the  target's  position,  (K  is  an  observed  target 

bearing  and  r^  is  the  range  of  the  initial  estimate  from  the 

station.  It  is  also  the  radius  of  a  circle  that  is  centered  on 

the  station  and  passes  through  the  initial  estimate  as  shown  in 

Figure  6  and  u^ ,  and  w^  are  related  arc  lengths  on  this  circle 

with  u.  =  w.  -V..  In  this  relation,  w.  is  known  and  v.  can  be 
111  1  1 

expressed  in  terms  of  the  unknown  coordinates  of  the  target  with 
an  approximation  that  does  not  involve  To  do  this,  consider 

a  rectangular  coordinate  system  whose  origin  is  at  the  position 
of  the  initial  estimate  and  whose  axes  are  oriented  like  those 
shown  in  Figure  6.  To  first  order  in  this  system,  v^  =  x  cos  3^^  - 
y  sin  where  x  and  y  are  the  unknown  target  coordinates  and 
-  4.  =  ( u .  -  c.)  -  (x  cos  3.  -  Y  sin  B.)/r..  The  use  of  this 

11  11  I"*  11 

relation  implies  that  the  bearing  line  determined  by  3^^ ,  the 

bearing  of  the  initial  estimate,  is  approximately  parallel  to  the 

true  bearing  line  determined  by  The  use  of  this  relation  for 

all  stations  imposes  the  constraint  on  the  maximum  likelihood  bearing 

estimates  by  replacing  what  would  otherwise  have  been  estimates 

of  n  independent  bearings  through  by  estimates  of  two  independent 

quantities,  the  rectangular  coordinates  x  and  y.  It  also  implies 


that  the  initial  estimate's  range  from  a  station  is  approximately 
the  target's  range  from  a  station,  that  is,  that  the  initial 
estimate's  position  is  relatively  close  to  the  target's  position. 


1 


Since  (0^  -  ())^)  ,  the  likelihood  of  observed  bearings 


0,  ,  0_,  0  can  be  written  as; 

1 '  2  n 


'  ®n)  = 


1  /Hr  a . 

1 


1  r  2  .2 

exp  -  j  E  u./o. 


where  a.  =  r.e.  and  e.  is  the  standard  deviation  of  0..  The 
1111  1 

maximum  likelihood  estimates  for  through  are  determined 

by  the  estimates  for  x  and  y  that  make  L(0j^,  ©2,  •••»  6^^)  a 

maximum.  In  this  case,  making  L(0j^,  02»  •••/  a  maximum  is 

^  2  2 

equivalent  to  making  I  (u./a.)  a  minimum.  So,  to  find  the  maximum 

1  ^  ^ 

likelihood  estimates  x  and  y,  solve  the  following  two  equations 

for  X  and  y: 


0  (  iln  L) 

3x 


=  0  and 


3(P,n  L) 
dy 


y  =  y 


X  =  X 


y  =  y 


With  the  constrai 


int  given  by  u .  =  w.  -  x  cos  0.  +  y  sin  0.  where 

^  -'ll  1  1 


w.  =  r.  (0.  -  0.),  the  equations  can  be  written  as  follows: 
1  11  1  ^ 


(w .  -  X  cos  6.  +  y  sin  3.)  (cos  3 . ) /o . 
1  1  ■*  1  11 


r  =  0 


n 

X  (w . 


-  X  cos  3.  +  y  sin  3.)  (sin  3.)/o  =  0, 


In  terms  of  the  following  quantities; 


E(cos^  6^)/o?  , 


E(sin^  S^)/a?  , 


E  =  E (w.  sin  6 ■ ) /of  , 
1  1  '  1 


the  equations  become: 


E(sin  8^  cos  8^)/o^  , 


E (w^  cos  , 


Ax  -  By  =  D 


Bx  -  Cy  =  E 


The  solutions  are: 


X  =  (BE  -  CD)/(B'^  -  AC) 


y  =  (AE  -  BD)/(B  -  AC) 


A  confidence  region  can  be  constructed  about  an  estimated 
position.  In  order  to  indicate  how  this  can  be  done,  a  proba¬ 
bility  region  about  the  true  position  will  be  considered  first. 

Both  X  and  y  are  values  of  random  variables.  If  a  new  set 
of  bearings  0^^,  02/  •••,  observed  (for  the  same  initial 

estimate  and  a  fixed  target),  in  general,  a  new  pair  of  values 
X  and  y  will  be  obtained. 


If  X  and  Y  represent  these  random  variables,  from  (1)  and  (2) , 

1  n  2 

(3)  X  =s  — i - E  (W./af)  (B  sin  -  C  cos  6.) 

(B'^-AC)  1  ^  ^  ^  ^ 

/si 

(4)  Y  -  — = -  E  (W./aT)(A  sin  B.  -  B  cos  B  • ) 

(B‘^-AC)  1  ^  ^ 

where  =  r^(0^  -  Bj^)  . 

Since  X  and  Y  are  a  linear  combination  of  the  n  normal  random 
variables  W2,  .../  W^,  or  equivalently  of  the  n  normal  random 
variables  0^^,  ©2,  0^,  they  have  a  joint  normal  distribution. 

Since  E(Wj^)  =  -  Bj^)  ,  if  B^^  =  for  i  =  1,  2,  n,  that 

is,  if  the  initial  estimate  of  the  target's  position  is  at  the 
target's  position,  E(Wj^)  =  0  for  i  =  1,  2,  ...,  n.  In  this  case 
E(X)  =0  and  E(Y)  =0  and  the  joint  normal  distribution  is  centered 
on  the  object's  position.  To  the  degree  of  the  approximations 
that  have  been  made,  this  is  also  true  if  the  initial  estimate  is 
not  at  the  target's  position. 

A  region  of  minimum  area  for  a  given  probability  of  contain¬ 
ment  of  an  estimated  position  can  be  determined.  The  region  is 
bounded  by  an  ellipse  which  is  centered  on  the  object's  position 
and  whose  axes  lie  along  the  axes  of  an  x'y' -coordinate  system 
that  is  obtained  by  rotating  the  xy-coordinate  system  that  is 
centered  on  the  object's  position  through  an  angle  y .  In  this 
system,  is  0,  that  is  X'  and  Y'  are  independent  normal  ran- 

dom  variables.  The  two  coordinate  systems  are  illustrated  in 


Figure  7.  The  coordinates  of  a  point  in  the  two  systems  are 
related  by 

x'  =  X  cos  Y  -  y  sin  Y 

y'  =  X  sin  Y  +  y  cos  y 

These  relations  imply: 

2  2  2  2  2 

(5)  cos  y  -  2a^^  cos  y  sin  y  +  sin  y  , 

X '  X  xy  y 

2  2  2  2  2 

(6)  sin  Y  +  cos  y  sin  y  +  cos  y 

y '  X  xy  y 

and 

2  2  2  2 

(7)  ^  =  (o^  -  0^)  sin  Y  cos  y  +  (cos  y  -  sin  y) 

X '  y  ’  X  y  xy 

where  y,  the  angle  of  rotation  of  the  coordinate  axes,  is  posi 

tive  in  the  clockwise  direction.  And  =  0  implies 

X  y 

tan  2y  = 

y  X 

With  the  initial  estimate  of  the  target's  position  at  the 

2 

target's  position  E(W^)  =  0  and  therefore  Var(W^)  =  for 

i  =  1,  2,  n.  In  this  case,  from  (3)  and  (4) 

^  1  2  2 
=  — 5 - 5"  ^  (1/c.)  (B  Sin  6.  -  C  cos  6.)^, 

X  (B^-AC)'^  1  ^  ^  ^ 


— > - *■  ^  (1/0.)  (A  sin  s.  -  B  cos  8..) 

(B'^-AC)'^  1  ^  ^  ^ 


1  "  2 

— X - E  (1/0.)  (B  sin  6.  -  C  cos  8^)  (A  sin  8-  -  B  cos  8-) 

(B^-AC)"^  1  ^  ^  ^  ^  ^ 


Using  the  definition  for  A,  B  and  C,  the  above  become 


^ 


X  (AC-B  ) 


^2  _  A 

^ ~  2  ' 
y  (AC-B  ) 


xy  (AC-B  ) 


So  tan  2Y  =  2B/(A-C)  for  8^^  =  <(>£f  i  =  1,  2,  n. 

With  the  target's  position  known  and,  consequently,  ((>^  known 

2  2 

for  i  =  1,  2,  ...,  n,  the  above  expressions  for  o^  ,  o^  ,  o  5.^  and  y 

X  y  xy 

can  be  used,  since  the  initial  estimate  of  the  target's  position 

can  be  taken  as  the  target's  position. 

With  values  for  0- i  o^  ,  o^^  and  Y,  values  for  o-,  and  o^, 

X  y  xy  X  y 

can  be  found  by  using  equations  (5)  and  (6).  The  probability 
that  an  estimated  position  will  be  within  an  ellipse  of  semiaxes 
ko^,  and  ko^ ,  which  is  centered  on  the  target's  position  is 


2 

1  -exp(-k  /2) .  This  result  can  be  found  by  integrating  the  bi¬ 
variate  normal  density  over  the  ellipse.  And  the  area  of  the 
2 

ellipse  IS  k 

Given  estimates  x  and  y  found  by  using  Equations  (1)  and  (2) , 

the  ellipse  with  semi-axes  ko'^,  and  ka^,  in  a  coordinate  system 

X  y 

that  is  centered  on  the  point  (x,y)  and  has  been  rotated  through 

2 

an  angle  y  is  a  1  -  exp(-k  /2)  confidence  region.  This  follows, 
since,  to  the  degree  of  the  approximations  involved,  the  bivar¬ 
iate  normal  distribution  of  X  and  Y  is  centered  on  the  target's 

position.  The  confidence  ellipse  is  defined  if  and  O''-' 

X  y  xy 

can  be  found,  that  is  if  the  elements  of  the  covariance  matrix 

can  be  found.  To  the  degree  of  the  approximations  involved, 

this  can  be  done  as  follows:  First,  assume  the  initial  estimate 

of  the  target's  position  is  at  the  target's  position.  Then, 

2  2 

values  for  ^  a*  ^  and  y  can  be  determined  by  using  Equations 

(8) ,  (9)  and  (10)  These  values  can  then  be  used  to  determine 
2  2 

,  0-',  and  by  using  Equations  (5),  (6)  and  (7).  Now, 

X  y  X  y 

with  a  value  for  k,  a  confidence  region  can  be  constructed.  To 
the  degree  of  the  approximations  involved,  the  shape  of  the  con¬ 
fidence  region  is  independent  of  both  the  target's  position  and 
of  the  initial  estimate  of  the  target's  position. 

For  the  case  where  bearings  are  taken  from  the  target  on  two 
or  more  stations,  0^^  is  the  reciprocal  of  the  bearing  taken  from 
the  target. 

A  discussion  of  the  theory  of  bearings  only  position 
estimation  procedures  for  situations  similar  to  the  one  considered 
here  is  given  in  Reference  2. 


The  following  equations  are  evaluated  in  the  program  to 
determine  the  coordinates  x*  and  y*  of  the  initial  estimate: 

X*  sin  (02-6]^)  =  C  Pj^  sin  (otj^-0j^)]  sin  02 

-  [  p2  sin  (QI2-02)]  sin  0^^ 

and 

y*  sin  (92“^^1^  ~  ^  ^1  (oij^-0j^)  cos  02 

-  L  P2  sin  {QI2-02)]  cos  0^^  . 

Using  the  point  determined  by  two  lines  of  bearing  as 
the  initial  estimate  was  suggested  by  a  similar  procedure 
described  in  Reference  3. 


PEST  is  a  progrcun  that  implements  a  position  estimation 
procedure  that  is  described  in  Appendix  1.  With  two  or  more 
observations  on  or  from  an  object  from  or  on  two  or  more 
stations,  the  program  will  generate  a  position  estimate  for  the 
object  with  an  associated  confidence  region.  After  program 
initiation,  the  program  user  is  presented  with  two  options;  (1) 
input  bearings  from  stations  on  an  object,  (2)  input  bearings  on 
stations  from  an  object.  To  choose  the  first  option,  press  T 
or  t.  To  choose  the  second  option,  press  S  or  s.  Bearing 
information  is  entered  in  the  following  order:  (1)  the  observed 
bearing  on  or  from  a  station,  (2)  the  bearing  of  the  station 
from  a  reference  location,  (3)  the  range  of  the  station  from  the 
reference  location  and  (4)  the  bearing  error.  The  reference 
location  can  be  a  station  location.  In  this  case,  the  station’s 
bearing  and  range  from  the  reference  location  are  both  taken  to 
be  zero.  After  two  observed  bearing  observations  and  their 
associated  station  and  bearing  error  information  have  been  input, 
the  user  is  given  the  option  of  (1)  generating  an  estimate  or 
(2)  continuing  to  enter  bearing  observation  data.  A  position 
estimate  for  an  object  is  specified  in  terms  of  its  range  and 
bearing  from  the  reference  location.  An  associated  elliptical 
confidence  region  which  is  centered  on  the  estimate  is  specified 
in  terms  of  the  length  and  direction  of  its  major  axis,  the 
length  of  the  its  minor  axis  and  the  probability  (confidence) 
that  it  contains  the  object.  The  program  user  is  given  the 
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option  of  defining  the  confidence  region  in  terms  of  (1)  size 
or  (2)  containment  probability.  For  example,  the  size  of  an 
ellipse  with  a  containment  probability  of  .8647  is  2.  In  this 
case,  the  ellipse  is  referred  to  as  a  two  sigma  ellipse,  since 
the  semi-major  and  semi-minor  axes  of  the  ellipse  are  two  standard 
deviation  units  in  length.  To  choose  the  first  option,  press 
S  or  s.  To  choose  the  second  option,  press  P  or  p.  The  estimate 
can  be  recalled  by  pressing  E  or  e.  Additional  observations  can 
be  entered  by  pressing  C  or  c.  To  quit,  press  Q  or  q. 

The  position  estimation  procedure  requires  an  initial 
estimate  of  the  object's  position.  In  the  program,  the  initial 
estimate  is  at  the  intersection  of  the  two  bearing  lines  that 
correspond  to  the  first  two  bearings  that  are  input  to  the 
program.  Because  of  this,  the  first  two  bearing  inputs  should 
be  from  or  on  the  two  stations  that  are  estimated  to  provide  the 
best  initial  estimate.  Note,  if  bearing  errors  are  large 
relative  to  the  angular  separation  of  the  two  stations  as 
seen  from  the  object,  the  bearing  lines  from  the  two  stations 
may  not  intersect.  If  this  is  the  case,  the  reciprocal  bearing 
lines  will  intersect  and  a  gross  error  in  the  final  position 
estimate  could  occur. 

The  procedure  which  is  essentially  equivalent  to  one 
described  in  Reference  2  is  based  on  the  following  assumptions: 

(1)  Bearings  are  taken  on  or  from  an  object  of  unknown  position 
from  or  on  two  or  more  stations  of  known  position.  (2)  The 
distances  involved  are  such  that  the  object  and  the  stations 


can  be  considered  to  be  located  on  a  plane  surface  (a  flat  earth) 
(3)  The  values  of  observed  bearings  on  or  from  stations  are 
determined  by  independent  normal  random  variables  each  with  a 
known  standard  deviation  (the  bearing  error)  and  a  mean  equal 
to  the  true  bearing  (zero  bias) .  (If  there  is  bias,  it  is  known 


and  removed . ) 


COMP  is  a  program  that  implements  a  procedure  for  combining 
a  group  of  position  estimates  for  an  object.  In  order  to  use  the 
procedure,  the  estimates  must  have  associated  elliptical 
confidence  regions.  When  used,  the  procedure  combines  the 
estimates  and  their  associated  confidence  regions  into  a  single 
composite  position  estimate  and  associated  elliptical  confidence 
region.  After  program  initiation,  the  program  user  is  prompted 
to  input  the  number  of  elliptical  areas  (confidence  regions)  to 
be  combined.  The  user  is  then  presented  with  two  options  for 
specifying  the  areas:  (1)  by  containment  probability  or  (2)  by 
size.  For  the  first  option,  the  areas  are  specified  by  the 
probability  (confidence)  that  they  contain  the  object.  This 
option  is  chosen  by  pressing  P  or  p.  For  the  second  option,  the 
areas  are  specified  in  terms  of  ellipse  size.  This  option  is 
chosen  by  pressing  K  or  Jc.  The  size  of  an  ellipse  with  a 
containment  probability  of  .8647  is  2 .  In  this  case,  the  ellipse 
is  referred  to  as  a  two  sigma  ellipse,  since  the  semi-major  and 
semi-minor  axes  of  the  ellipse  are  two  standard  deviation  units 
in  length.  The  orientation  of  an  ellipse  is  specified  in  terms 
of  the  direction  (the  angle  delta  in  the  program)  of  its  major 
axis.  The  position  estimates,  the  centers  of  the  elliptical 
areas,  are  specified  in  terms  of  latitude  and  longitude.  The 
user  is  given  the  option  of  determining  a  confidence  region  in 
terms  of  (1)  its  containment  probability  (confidence)  or  (2)  its 
size.  To  choose  the  first  option,  press  P  or  p.  To  choose  the 


second  option,  press  K  or  k.  To  recall  the  position  estimate, 
press  E  or  e.  To  quit,  press  Q  or  q. 

The  procedure  which  is  described  in  detail  in  Section  III  is 
based  on  the  following  assumptions!  (1)  The  individual  estimates 
are  determined  by  independent  normally  distributed  random 
vectors  with  known  covariance  matrices  and  a  common  but  unknown 
mean  vector  whose  components  are  the  object's  rectangular 
coordinates.  (2)  The  object  is  located  in  the  plane  of  the 
coordinate  axes,  a  plane  tangent  to  a  spherical  earth.  In 
the  program,  the  coordinates  of  the  first  entry  are  the  coordi¬ 
nates  of  the  point  of  tangency.  Because  of  other  uncertainties, 
this  flat  search  assumption  should  not  introduce  significant 
estimation  errors  for  the  distance  scale  for  which  the  program 
is  intended.  In  a  sense,  the  procedure  is  a  generalization  of 
the  procedure  described  in  Reference  4  which  is  limited  to 
circular  confidence  regions. 


10  CLSi  SCREEN  0,0:  DIM  A (7) 

20  PRINT  "T^bsarings  on  the  target,  “i  PRINT  ''S"bearings  on  the 
stations. " 

25  A<=INKEY* 

30  IF  At="T"  DR  A»»"t"  THEN  JJ-Oi  GOTO  60 

40  IF  A*=''S"  OR  A*»''s"  THEN  JJ-ll  GOTO  60  ELSE  GOTO  25 

60  CLS: INPUT  "observed  bearing";Pt  P»P*ATN< 1 ) /45;  IF  JJ»1  THEN 

p=p+iao 

70  INPUT  "station  bearing"|Q:  Q»Q*ATN < 1 > /45: INPUT  "station 
range";R:  INPUT  "bearing  error";0:  0»0*ATN( 1 ) /45 
80  IF  I«2  THEN  GOTO  140 

90  1=1+1:  A(I-1>=P:  A(I+1)=Q:  A(I+3)»R:  A(l+5)»0:  IF  1=1  THEN 
GOTO  60 

100  X=A(4)*SIN(A(2)-A(0> ) :  Y=A (5) *SIN ( A <3> -A < 1 > > :  Z=SIN ( A ( 1 > -A (0) 
):  IF  Z=0  THEN  GOTO  330 

no  U=(X*SIN(  A(1))-Y*SIN  (A(0)))/Z:  V=(X*COS  ( A ( 1 > > -Y*COS (  A(0>> 
)  /Z 

120  FOR  M=0  TO  1 

130  P=A(M):  Q=A(M+2>:  R=A(M+4>:  0=A(M+6>:  GOSUB  390:  NEXT  M:  GOTO 
150 

140  GOSUB  390 

150  CLS:  PRINT:  PRINT  "E=Est  C=Cont" 

160  A»=1NKEY4:  IF  A4="E"  OR  A*="e"  THEN  GOTO  IBO 
170  IF  A4="C"  OR  A4="c"  THEN  GOTO  60  ELSE  GOTO  160 
ISO  CLS;  F=(B*B-A*C):  IF  F=0  THEN  GOTO  380 

190  X1=U+(B*E-C*D) /F:  Y1=V+ ( A*E-B*D) /F:  GOSUB  440:  K-Rl:  J»B1 
200  T=SGN(B)*ATN(1) :  IF  A=C  GOTO  220 
210  T=.5*ATN(2*B/ (A-C) ) 

220  G=<C*C0S<T)*C0S<T)-2*B#C0S<T>*S1N<T)+A*SIN(T)#S1N(T>  >/-F: 
G=SQR(G) 

230  H=  <C*SIN  <T) *SIN (T) +2*B*C0S (T) *S1N  <T) +A*COS (T>  *COS (T>  > /-F: 

H=SQR(H):  IF  H>=G  GOTO  250 

240  Z=H:  H=G:  G=Z:  T=T+2*ATN(1) 

250  CLS:  PRINT;  PRINT  USING  \####. ##" 5 "bearing=" , J*45/ 

ATN ( 1 ) 

255  I1=K;  GOSUB  500;  PRINT  "range»"jll:  PRINT 
260  PRINT  "S=Size  P=Prob  E«Est  C=Cont  Q-quit" 

265  A»= INKEY* 

270  IF  A*="S"  OR  A*="s"  THEN  CLS:  GOTO  320 

280  IF  A*="F"  OR  A*="p"  THEN  CLS:  GOTO  350 

290  IF  A*="E"  OR  A*="e"  THEN  CLS:  GOTO  250 

300  IF  A»="C"  OR  A*="c"  THEN  CLS:  GOTO  60 

310  IF  A*="Q"  OR  A*="q"  THEN  END  ELSE  GOTO  265 
320  INPUT  "size";S:  IF  S<=0  THEN  GOTO  320 
330  0=l-EXP(-S*S/2) 

340  PRINT  USING  "\  j "probabi 1 i ty«" ,0:  GOTO 

370 

350  INPUT  "probability";0:  IF  0>=1  OR  0<=0  THEN  GOTO  350 
360  S*SQR(-2*L0G<l-0) ) :  PRINT  USING  "N  ##" i "si ze«" ,S 

370  X=S«G:  I1=2*X:  GOSUB  500:PRINT  "major  axis*"; II 

371  N=T*45/ATN(1) :  IF  N<0  THEN  N=N+180 

372  PRINT  USING  "\  \####.«#"; "direction*" ,N 

375  Y«S*H: I1«2*Y:  GOSUB  500:  PRINT  "minor  axis-"; II;  I1=4*ATN(1)* 
X*Y:  GOSUB  500:  PRINT  "area-"; II:  GOTO  260 
380  PRINT  "no  solution":  END 


390  X1»IJ-R*SIN(Q)  I  Yl«V-R«C0S(i>i«  I  BOsUB  440 
400  W-P  Bli  L»Rl*Oi  IF  L-0  THEN  GOTO  380 

410  G-CI)S(B1)/Li  H«SIN(B1)/Li  IF  W>-4«ATN<1>  THEN  W-W-8*ATN ( 1 )  i 
GOTO  43<  > 

420  IF  i4<— 4*ATN(1>  THEN  W«W <-8*ATN ( 1 ) 

430  W-W/O:  A»6«G-i-At  Ba6*H-^Bt  C-H«H-i-C:  D->W«G-i-Di  E-W«H-«-Ei  RETURN 
440  R1»3QR(X1*X1+Y1*Y1) :  IF  Rl-0  THEN  Bl»Oi  RETURN 
450  IF  f^BS(Xl/Rl>=l  THEN  Ml-SGN<X1 )*ATN(1)*2  ELSE  M1»ATN<X1/R1/ 
SQR(1-X 1*X1/R1/R1> ) 

460  IF  i^BS(Yl/Rl)=l  THEN  BI=2*ATN ( 1 ) * ( 1-SGN<Y1 ) )  ELSE  Bl=2*ATN(l) 
-ATN(Yl/Rl/SaR(l-Yl*Yl/Rl/Rl ) ) 

470  IF  I11<0  THEN  B1=8*ATN  ( 1 ) -B1 
480  RETURN 

500  11-100*11:  Il-INT<Il)t  Il-Il/lOO 
510  RETURN 


COMP 


5  CLS:  PI«4*ATN(l):  INPUT  "Numbar  of  Elliptical  Ar««a'‘|AOi  PRINT 

10  PRINT  "Ar*a  Definitions! "iPRINT 

11  PRINT  "By  Containment  Probability,  Press  P. PRINT 

12  PRINT  "By  Sigma  Size,  Press  K." 

15  E»= INKEY* 

16  IF  E*="P"  OR  E**"p"  OR  E*="K"  OR  E*»"k"  THEN  GOTO  17  ELSE  GOTO 
15 

17  CLSi  PRINT  "For  the  Latitude  and  Longitude  entry  format,  press 
F.  OtherMise,  press  C. " 

IS  K*=INKEY* 

19  IF  K*="C"  OR  K*»"c"  THEN  GOTO  26 

20  IF  K*="F"  OR  K*»"f"  THEN  GOTO  23  ELSE  GOTO  IQ 

23  CLS:  PRIN1  "Latitude  and  Longitude  are  entered  in  degrees  and 
minutes  and  tenths  of  minutes  in  the  form:  DDD-MM.liX  where  X  is 
N,S,W  or  E.  Leading  zeros  are  optional,  but  the  -  must  be 
included.  Press  C  to  continue." 

24  C*= INKEY* 

25  IF  C*="C"  OR  C*="c"  THEN  GOTO  26  ELSE  GOTO  24 

26  CLS! FOR  0=1  TO  AO 

27  EE=0:  PRINT 

28  INPUT  "LAT"5A*!  GOSUB  410:  IF  EE=1  THEN  GOTO  27 

29  Y=VA 

30  INPUT  "LONG"; A*:  GOSUB  410i  IF  EE=1  THEN  GOTO  27 

31  X=VA 

32  IF  0>1  THEN  GOTO  34 

33  Y0=Y!  X0=X!  C0=C0S(Y0*PI/1Q0) 

34  Y=(Y-Y0)*60:  X=(X0-X) 

35  IF  X>270  THEN  X=X-360 

36  IF  X<-270  THEN  X=X+360 

37  X=X*60*C0 

38  INPUT  "Delta";S 

40  IF  S>=180  OR  S<0  THEN  GOTO  35 

45  INPUT  "Major  Axis"  ;P:  INPUT  "Minor  Axis";Q!  P-P/2!  Q-Q/2 
50  IF  E*="K"  OR  E*="k"  THEN  GOTO  60 
55  INPUT  "P";P3!  R=SQR<-2*L0G(1-P3) ) !  GOTO  70 
60  INPUT  "K"jR 

70  IF  S>=135  THEN  S=S-180iT»P  !P=Q!  Q-T:  GOTO  110 

80  IF  S>=45  THEN  S=S-90!  GOTO  110 

90  IF  S>=0  THEN  T=P!  P=Q!  Q=T!  GOTO  110 

110  P=P/Ri  Q=Q/R 

120  S=S*PI/180 

1 30  U=P*P*COS ( S ) *COS ( S ) +Q*Q*S I N  <  S ) *S I N  <  S ) 

1 40  V:=P«P*S  I N  <  S  >  «S  I N  ( S )  -*-a*Q*COS  ( S  >  *COS  <  S ) 

150  W=(Q*Q-P*P)*SIN(S)*COS(S) 

160  Z=U*V-W*W:  A-A+V/Z!  B-B-W/Zi  C-C+U/Zi  D=D+V/Z*X-W/Z*Y 

170  E=E-W/Z*X>U/Z*Y!  F-F+V/Z*V/Z*U!  G=G+W/Z*W/Z*U 

180  H=H-V/Z*W/Z*U!  1=1+W/Z*W/2*V:  J»J+U/Z*U/Z*V 

190  K=K+W/Z*U/Z*V!  L=L-V/Z*W/Z*W 

200  M=M+<V/Z*U/Z+W/Z*W/Z)*W!  N=N-W/Z*U/Z*W 

210  NEXT  0 

215  CLS 

220  0=A*C-B*B!  X*<C*D-B*E) /O:  Y* < A*E-B*D> /O 

230  R=(C*C*(F+I+2*L)-2*C*B*(H+K«-M)+B*B*<G+J-«-2*N)  >  /  (0*0) 
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240  S*  <  A*A*  <  G+ J+2*N ) -2*A*B# ( H+K+M ) +B*B* ( F+ 1 +2*L ) ) / <  0*0 ) 

250  T»  <  ( A*C-«-B*B  >  *  ( H+K+M )  -C*B*  ( F+ 1 +2*L  )  -B*A*  ( G+ J +2*N )  )  /  <  0*0 ) 

255  Y-YO+Y/60:  VA-Yi  GOSUB  470i  GOSUB  500 
257  PRINT  “LAT  -  "jL* 

25B  XaX0-X/(60*C0) I  VA»Xi  GOSUB  470i  GOSUB  540 
259  PRINT  "LONG  »  '’|G* 

270  A«SGN<T)*ATN(1) :  IF  S«R  GOTO  280 
275  A*.5*ATN(2*T/ (S-R) ) 

280  B«R*COS ( A ) *COS ( A > -2*T*C0S ( A ) *S I N ( A >  +S*S IN(A)*8IN(A) 

290  C=R*SIN(A)*SIN(A)+2*T*SIN<A)*CQS(A)+S*C0B(A>*C0S(A) I  A»A*180/ 
PI 

300  IF  B>=C  THEN  D=A+90t  GOTO  350 

305  E=B:  B=C!  C=E:  IF  A<0  THEN  D«A+180t  GOTO  350 

310  D*A:  GOTO  350 

320  PRINT:  INPUT  "P";P3:  K=SQR (-2*L0G < 1-P3) ) i  Fl»li  GOTO  330 
325  PRINT:  INPUT  "K";K:  P3=1-EXP (-K*K/2) :  F1»0 
330  CLS:  PRINT  "Major  Axis  =  ";2*K*SQR(R) 

335  PRINT  "Minor  Axis  =  ";2*K*SQR(S) 

340  PRINT  "Delta  =  ";D:  IF  Fl=l  THEN  GOTO  345 

341  PRINT  "P  =  ";P3:  GOTO  350 
345  PRINT  "K  =  ";K 

350  PRINT:  PRINT  "E  =  Est  P>=P  K»K  Q»  Quit" 

355  A4=INKEY4 

360  IF  A4="E"  OR  A4="e"  THEN  GOTO  390 

365  :F  A4="P"  or  A*="p"  THEN  GOTO  320 

370  IF  A*="K"  OR  A«="k"  THEN  GOTO  325 

380  IF  A4*"Q"  OR  A4= "q"  THEN  GOTO  400  ELSE  GOTO  355 

390  CLS:  PRINT  "LAT  =  ";L»:  PRINT  "LONG  -  "jG4:  GOTO  350 

400  CLS:  END 

410  VA=VAL<A*) 

415  D4«MID4<A4,4,1) -  IF  D4="-"  THEN  GOTO  435 

420  D4=MID*<A4,3,1) :  IF  D»="-"  THEN  GOTO  430 

423  D»=MID4<A4,2,1) :  IF  D4="-"  THEN  GOTO  425 

424  CLS:  PRINT  "Data  entry  error,  restart  entry  for  the  area.": 
EE=1:  RETURN 

425  M4=MID4(A»,3) :  GOTO  440 

430  M»=MID4(A*,4) :  GOTO  440 

435  M*=MID»(A*,5) :  GOTO  440 

440  VM=VAL<M4):  VA=VA+VM/60 
450  R*=RIGHT*(M*,1) 

460  IF  R*="S"  OR  R4*"s"  OR  R**"E"  OR  R4="e"  THEN  VA*-VA 
465  RETURN 

470  QO=SGN(VA):  VA»ABS<VA) 

473  IF  VA>180  THEN  VA=*(VA-360)  ELSE  GOTO  480 
475  QQ=-QQ:  VA=ABS<VA) 

480  DD*FIX(VA):  FF^VA-DD:  MM«-F I X < 6000*FF > / 1 00 
490  RETURN 

500  IF  QQ=1  THEN  Q4»"N":  GOTO  520 
510  IF  QQ»-1  THEN  Q4*"S"  ELSE  Q4-"  " 

520  L»=STR»<DD)+STR4<MM)+Q4 
530  RETURN 

540  IF  QQ=1  THEN  Q**"W":  GOTO  560 
550  IF  QQ»-1  THEN  Q4»"E"  ELSE  Q»-"  " 

560  G#=STR4<DD)+STR4<MM)+Q* 

570  RETURN 
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